Michelson interferometry with quantum noise reduction 
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A Michelson interferometer with noise reduction to sub-shot noise levels is proposed and realized. 
Multiple measurements of a single signal beam are taken and the quantum property of light plays an 
essential role in the principle underlying this interferometry. The method makes use of the coherent 
state of light and requires only a simple modification to the standard Michelson interferometer. 
The surface fluctuation spectra of liquids are measured using this method down to a few orders of 
magnitude below the shot noise level. The spectrum derived from hydrodynamical considerations 
■ agrees well with the observed results for water. However, for oil, slight deviations are seen at high 

frequencies (> 1MHz), perhaps indicating its more complex underlying physics. The measurement 
requires a relatively low light power and a short time, so that it has a wide range of applicability. 
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Interferometry is a most precise tool for detecting small displacements and hence is used in a very broad range 
of areas in physics, from microscopic spectroscopy such as ours, structural measurements of optical elements to 
astrophysics [lj and arguably the most sensitive measurement, the attempts to detect gravitational waves 
Surface thermal fluctuations have been measured with interferometry previously, on mirrors using high power 
lasers [|| and on liquid drop surfaces attached to fiber tips [Hi Spectral properties of surface waves on simple 
liquids have also been measured by using them as gratings |g, 0] . However this approach is difficult to apply 
to dissipative liquids such as oil, since they do not create well defined waves that act effectively as gratings. 
Spectra of surface inclination fluctuations have been measured (54ll| using the surface as an optical lever. Also, 
exceptionally large surface fluctuations due to low surface tension have been observed using other methods 

In this work, we measure surface thermal fluctuation spectra of simple liquids over a wide frequency range 
(1kHz ~ 40 MHz) and down to few orders of magnitude below the shot noise level, using Michelson interfer- 
ometry. While these spectra reflect fundamental physics principles of liquids and are interesting from a physics 
perspective^ 0, EH, the spectra we obtain have not been done so previously, to our knowledge. The main reason 
for this is that these fluctuations are small and are buried under the shot noise. It is usually believed that the 
i-C ' shot noise sets the limit for the signal-to- noise ratio that can be obtained under normal circumstances [21 15. [llj . 

with the exception being measurements involving sub-Poissonian photon statistics. We show both theoretically 
and experimentally that the standard Michelson interferometer, when combined with the quantum nature of 
light can achieve signal-to-noise ratio unlimited by the shot noise. In essence, while both thermal fluctuations 
and shot noise are random, the former are classical in their origin and the latter has fundamentally quantum 
nature and their difference enables us to separate these two. By being able to analyze the thermal fluctuation 
spectra to such precision, their theoretical understanding can be examined in detail. The traditional hydrody- 
namical description of simple liquids works well for water, but some deviations from it for liquids such as oil 
are observed at high frequencies. 

We now summarize briefly the general principle behind the noise reduction employed here, which is not 
limited to optical measurements nor to shot noise [l(| [H|: In analyzing weak signals, if the signal has a definite 
periodicity, we can accumulate data in accordance with the period to suppress the noise. However, in cases 
where there is no definite periodicity, it is difficult to separate out random signals, such as thermal fluctuations, 
from random noise. In fact, given a single measurement D\ = S + N± of a random signal, there is no way, 
even in principle, to distinguish the signal S from the noise Ni, which inevitably occurs. Shot noise is a typical 
example of such random noise. To overcome this obstacle, we perform an additional measurement of the same 
signal Z?2 = S + N2 at the same time, whose noise N2 is independent of Ni . Then, by taking the correlation of 
D1.2 and averaging over time, we obtain 

(iT^) ^ (\S\ 2 ) (AA^oo) , (1) 

where J\f is the number of averagings and tildes denote Fourier transforms [Tol fll|. While the above principle 
is simple, the crucial point is to arrange multiple measurements of the same signal in such a way as to ensure 
that the noise in them are uncorrelated. 

To measure surface fluctuations, laser light is shone on the surface and the reflected light is used as a signal 
for measurement. The quantum property of light plays a critical role in our experiment; the beam splitter 
randomly partitions the photons from the single signal beam to one of the two detection systems (Fig. HJa)). 
The photons giving rise to the shot noise is random, when a coherent light source is used. Therefore, the 
noise in the two detectors are uncorrelated and is eliminated when we compute their correlation, Eq. ([1]). It 
is important to note that splitting the beam by itself does not guarantee the independence of the noise in the 
two measurements. Had we used a squeezed light source, for instance, the photons in the two measurements 
would have been correlated, so that the noise could not have been eliminated by using the correlation of detector 
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measurements. Therefore, our approach is in contrast to those that use squeezed light sources to obtain sub-shot 
noise measurements 0, [H, EH- Using squeezed light sources, a reduction of the shot noise by a factor of two 
have been achieved and a light source squeezing factor close to 20 has been attained 17 1 . 
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FIG. 1: (a) By taking the correlation of two photodetector (PD1,2) measurements, we separate the signal from the 
uncorrelated random noise. (b)To eliminate AM noise that inevitably occur, we additionally use the differential detection 
of random signal 1 and 2, by adding PD3,4. Random signals 1,2 can either be phase inverted or uncorrelated. 



The signals obtained from the photodetectors (Fig. HJa)), in practice, additionally contain amplitude mod- 
ulation (AM) noise from the light source. This is common to the two measurements and hence will not be 
eliminated by taking their correlation. To reduce this noise, we further employ differential detection (Fig.[TJb)). 
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FIG. 2: Experimental setup that incorporates noise reduction due to correlation and differential detection shown in 
Fig. [1] NPBS: nonpolarizing beam splitter. 



Our Michelson interferometry setup that incorporates these principles is shown in Fig. [2] The light source for 
this experiment is a solid-state diode-pump laser (SAPPHIRE SF-532-50, Coherent) with a wavelength of A=532 
nm arranged so that the power is 500 fiW each at the sample and at the reference mirror. A typical measurement 
we show below requires a few seconds of measurement time. To implement the differential detection, a Faraday 
rotator is used for separating out the light reflected back from the sample at the polarizing beam splitter (PBS). 
The lengths of the signal and reference arms in the Michelson interferometer are adjusted to be equal within 
1 /mi, since a larger difference increases frequency modulation noise effects. To obtain maximal sensitivity in 
detecting sample surface displacements, the reference arm length is adjusted dynamically by a piezo-electric 
transducer (PZT) with feedback. In the Michelson interferometer, the light at the reference mirror and that at 
the sample need to have the same properties or the destructive interference within the beam spot will reduce 
the signal significantly. Light reflected by the reference mirror and the sample is detected by four photodiodes 
(S5973 Hamamatsu Photonics, Japan). The signal currents obtained through differential detection are digitized 
by an analog to digital converter (ADC, 8 bit, 125 Ms/s, PicoScope 5203, Pico Technology). Their Fourier 
transforms and the correlation in Eq. ([1} are calculated using a computer. 

Two light sources and two sets of detectors can also be used to obtain two independent measurements and 
extract signals at sub-shot noise levels from their correlation. This was the approach used for inclination fluctu- 
ation spectra of surfaces [13, HU- I n this work, by using a single light source, we make minimal modifications to 
the classic Michelson interferometer, essentially by just adding an additional beam splitter at the photodetector, 
to attain sub-shot noise measurements. This leads to a simple elegant setup. Using a single light source is not 
only simpler, but has an important practical advantage; to achieve maximal sensitivity in the interferometer, 
the path lengths need to be adjusted according to the wavelength of the light source, which is difficult to attain 
for multiple wavelengths simultaneously. 

The spectral function of thermal surface fluctuations for a simple liquid is determined from hydrodynamical 
considerations by its density p, surface tension a and viscosity r\ and is 

(2) 

where u = puj/{2rjk 2 ), y = pa/(4r] 2 k)^U^. Here, k,uj are the wave number and the angular frequency of the 
surface wave. Through Michelson interferometry, we measure fluctuations perpendicular to the liquid surface 
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and its spectrum can be computed analogously to the inclination fluctuation spectrajlOl lll| as 

poo 

5 h (/) = 2 / dkke- b2k2 / s P(k,2irf) , (3) 
Jo 

Here, / = lo /(2tt) and b is the beam diameter. Surface waves with wavelengths larger than the sample size 
are cut off so that the effects of gravity can be ignored. Also, surface waves with wavelengths smaller than the 
beam spot size is suppressed due to averaging. In the measurements below, the size of the surface sample is 
2.2 mm in diameter and b = 0.96 /Ltm. 
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FIG. 3: (Color online) Experimentally observed surface height fluctuation spectra for water (red) and ethanol (blue). 
Water has smaller fluctuations. Respective theoretical spectra are also shown (black, dashed), which agree with the 
experimental results and are almost invisible. For comparison, observed data for a single differential detection without 
using the correlation Eq. flj is shown for water (green) which is clearly dominated by the shot noise at higher frequencies. 
The theoretical value for the shot noise level, Eq. Q, is also indicated (black, dotted). 
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FIG. 4: (Color online) Experimentally observed surface fluctuation spectra for oil at various temperatures (red, green, 
blue, magenta from lower to higher temperatures). Fluctuations are larger at higher temperatures. Their respective 
theoretical spectra, Eq. ([3]), are also shown (black, dashed). (Inset)Temperature dependence of the deviation from 
theory of the fall off at high frequencies (> 1 MHz). 

We compare the measured spectra against their theoretical predictions, Eq. ([3]), for water and ethanol 
in Fig. [3J The physical properties of water, ethanol are well known for a given temperature and are 
(/9[kg/m 3 ],CT[kg/s 2 ],r/[kg/(m-s)],T[K]) = (1.0 x 10 3 ,7.5x 10~ 2 ,1.7 x 10- 3 ,275), (0.80 x 10 3 ,2.3 x 10- 2 ,1.7 x 
10 _3 ,278), respectively. Given the temperature, these properties and the spectral function P(k,u>) are com- 
pletely determined. The integrated spectrum Sh(f), Eq. (j3|), is uniquely determined by those properties and 
the beam size. Water and ethanol are cooled to suppress evaporation which can cause the sample surface to 
leave the focal plane. For water surface fluctuations, the theoretical formula derived from hydrodynamic con- 
siderations agree essentially perfectly with the observed spectrum. These measurements confirm the theoretical 
spectrum over a wide frequency range. For ethanol, good agreement is seen, except for the slightly slower fall 
off at frequencies above 10 MHz. We come back to this issue for the case of oil surface fluctuations. 
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e ~^r^2~7 — A/ , (4) 



Oil (Olympus immersion oil AX9602) surface thermal fluctuation spectra were obtained at various tempera- 
tures and compared to the theoretical spectra Eq. ^ in Fig. 01 There is a qualitative difference from the spectra 
of water and ethanol surface fluctuation spectra in Fig. [31 due to the high viscosity of the liquid. In particular, the 
spectral density decays much more slowly than that of water and oil at higher frequencies. For the oil, the tem- 
perature dependence of a, rj are not known so that they had to be deduced from the spectra and were found to be 
(p [kg/m 3 ],cr[kg/s 2 ],?7[kg/(m- s)],T[K]) = (0.92 x 10 3 , 3.2 x 10" 2 , 0.15, 294), (0.91 x 10 3 ,3.1 x lO" 2 , 0.096, 303), 
(0.89 x 10 3 ,2.9 x 10" 2 , 0.037, 332), (0.88 x 10 3 ,2.7 x 10" 2 , 0.025, 347). a, rj are determined to ten percent from 
the spectra. For oil, we see that the fluctuation spectrum changes with the temperature, as it should. Most of 
this change is caused by the decrease in the viscosity with higher temperatures, making the fluctuations larger. 

While the fluctuation spectra of liquids we study are usually regarded as well understood, the observed 
higher frequency (> 1 MHz) fall offs in the spectra for oil and ethanol are slightly slower than their theoretical 
predictions, especially at lower temperatures. This is possibly due to the more complex nature of the liquids 
that can not be explained just from the hydrodynamical considerations of simple liquids. The slower fall off 
of the spectrum has also been observed for solid materials and a gradual transition to such a dependence was 
seen for complex fluids 10]. We obtain the dependence Sh(/) ~ f a for higher frequencies (10 6 — 10 7 Hz) and its 
deviation 8a from the theoretical value in Eq. ([3]). In Fig.|4jinset), the temperature dependence of 8a is shown. 
The dependence can be reasonably well described by 8a = C exp([7 /k B T) with U = 7kJ/mol. This energy scale 
U is comparable to the latent heat for oil, which is consistent with more complex behavior, such as molecular 
interactions causing visible effects in the spectrum. 

In Fig. [3l we also included the results from one differential measurement of water surface fluctuations, in 
which the shot noise level is clearly visible. Theoretically, the shot noise level in the spectrum is 

A 2 1 + r 2 
32^ rIp D 

when the sensitivity in the interferometer is maximal. Here, r is the ratio of the reflectivities of the sample 
surface and the reference mirror, Ipd the signal photocurrent and e the electron charge. As can be seen in 
Fig. [3] and 01 shot noise is clearly eliminated when the correlation of the two differential measurements are 
taken. The observed shot noise level in Fig. [3] is roughly twice its theoretical value, Eq. This indicates that 
the interferometer sensitivity is not at its theoretical maximum, whose most likely cause is the aberration of 
the objective lens. 

Let us briefly describe the physical properties of the spectra Eq. (|3j: When IQ^/^-Kirf f / (pa 2 ) > 1, the liquid 
can be regarded as being highly viscous so that any liquid is dissipative at high enough frequencies. In the 
frequency region we study, water, ethanol have low and oil has high viscosity. The spectra have qualitatively 
different / dependence for low and high viscosities, as explained belowflollTTl]. In both cases, Sh(/) ~ k B T/(af) 
at low frequencies. For a liquid with low viscosity, the spectrum crosses over at / ~ \J a / (pb 3 ) to Sh.(f) ~ 
k B Ti]/ '(p 2 b 5 / 4 ). For a liquid with high viscosity, the spectrum crosses over at / ~ cr/ (rib) to Sh(/) ~ k B T / (rjbf 2 ). 
While outside the region of our measurements, at even higher frequencies, / > rj/(pb 2 ), the spectrum changes 
to Sh(f) ~ k B Tr]/(p 2 b 5 f ). For both low and high viscosities, the behavior of Sh(/) is dominated by a at 
low frequencies and is independent of rj since the dynamical time scale is relatively large. At high frequencies, 
the behavior is governed by rj and is independent of cr, to leading order. These crossovers seen in the spectra 
are similar to those in the surface inclination fluctuation spectra and arise from the properties of P(k, 27rf), 
Eq. © [10, [H| ■ However, the spectra differ qualitatively. In particular, Sh(f) has a distinctive 1// dependence 
at lower frequencies and is independent of the beam diameter b to leading order. Consequently, the height 
fluctuation measurements are more sensitive than the inclination fluctuation measurements at lower frequencies. 
This is natural since the inclination fluctuations due to longer wavelength fluctuations are smaller for a fixed 
beam size. At high frequencies, the cutoff frequency determined by b becomes important, since it is responsible 
for suppressing the shorter wavelength modes. 

In this work, we have proposed and implemented the Michelson interferometer with noise reduction at sub-shot 
noise levels. We make multiple measurements of a single signal light. The shot noise in them are independent 
due to their quantum nature, hence the noise can be eliminated through their correlation and the signal below 
this noise level can be extracted. While both random, the distinct difference between the thermal fluctuations 
and the shot noise, which are classical and quantum in their origins, allow for their separation. It would 
be interesting to understand the behavior when the method is applied to quantum fluctuation measurements 
with shot noise. At present, some applications such as quantum cryptography and quantum random number 
generation exist, yet there have not been many cases where the quantum property of light has been of practical 
value flij and it is satisfying to find it provide a practical advantage here. The measurement requires only a 
small light power, a relatively short time and a small sample surface (500 ^W, few seconds, diameter 1 /zm, in 
this work) . Our noise reduction can also be added to higher power measurements such as those in Q to reduce 
the noise further. Longer measurement times can lead to higher resolution, more averaging and less statistical 
error. Intcrfcromctry is perhaps the most commonly used method for precision measurements and shot noise 
can often be the major limiting factor in its accuracy 0, 0, [lij]. We believe that our method can be applied in 
various situations in which precision, low power or short measurement time is required. 
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Let us discuss the limitations of our approach: The approach is effective when the number of averagings 
Af = A/T is large, where A/ is frequency resolution in the spectrum and T is the total measurement time. The 
relative error in the spectrum is 1/\M. Af needs not be the same across the spectrum and if a same relative 
resolution is used, this statistical error is smaller at higher frequencies. While this integration is effective for 
stationary signals, the situation is more subtle for transient signals such as gravitational waves radiated from a 
single event. In such situations, the restriction that the measurement time needs to be within the time when 
the signal is present can be demanding. The duration of the signal needs to be long enough compared to the 
inverse of the desired frequency resolution, to reduce the noise using the correlation of measurements. 
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